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Abstract* 



A new model for the postcollapse equilibrium structure 
of virialized objects which condense out of the cosmo- 
logical background universe is described and compared 
with observations and simulations of cosmological ha- 
los from dwarf galaxies to X-ray clusters. The model 
is based upon the assumption that virialized halos are 
isothermal, which leads to a prediction of a unique non- 
singular isothermal sphere for the equilibrium structure, 
with a core radius which is approximately 1/30 times the 
size and a core density which is proportional to the mean 
background density at the epoch of collapse. These pre- 
dicted nonsingular isothermal spheres are in good agree- 
ment with the observations of the internal structure of 
dark-matter-dominated halos from dwarf galaxies to X- 
ray clusters. 



* Poster paper presented at the Institute for Theoretical 
Physics Conference on Galaxy Formation and Evolution, 
March 14-17, 2000, at the University of California at 
Santa Barbara. 



Equilibrium structure of virialized Inalos 



[Shapiro, Iliev, and Raga 1999, MNRAS, 307, 203 (Einstein- 
de Sitter case); Iliev and Shapiro 2000 (low-density, 
open and flat with cosmological constant)] 



• Problem and Motivation 

— Question: What equilibrium structure forms when 
a density perturbation collapses out of the ex- 
panding background universe and virializes? 

— An analytical model for the structure (e.g. mass 
profile, temperature, velocity dispersion, radius) 
of virialized cosmological halos would be a valu- 
able tool for the semi-analytical modeling of 
galaxy and cluster formation in a hierarchical 
clustering model like CDM. 

* Earlier work adopted crude approximations in- 
volving either uniform spheres or singular isother- 
mal spheres which resulted from top-hat per- 
turbation collapse 

* What is a more realistic outcome, even for 
the simple top-hat problem? 

— N-body simulations of CDM predict dark matter 
halo profiles with singular density profiles, but a 
finite density core is required to explain: 

* Dwarf galaxy rotation curves 

* Cluster mass profiles inferred from gravita- 
tional lensing 

— As a result, the cold, collisionless nature of CDM 
has recently been re-examined to allow for vari- 
ations which affect the post-collapse equilibrium 
structure of halos. 



— Suppose we ignore the details of this relaxation 
process and adopt the assumption that the final 
equilibrium is isothermal. 

— Solve this problem and compare the result with 
dwarf galaxy rotation curves and X-ray cluster 
data. 

• Model: 

— Top-hat density perturbation collapses and viri- 
alizes 

— Virialization leads to a truncated isothermal sphere 
in hydrostatic equilibrium (TIS) solution of 
the Lane-Emden equation (modified for A 7^= 0) 

— Total energy of top-hat is conserved thru col- 
lapse and virialization 

— Postcollapse temperature set by virial theorem 
(including effect of finite boundary pressure) 

• Is the solution uniquely determined? - No, some 
additional information is required: 

1) Minimum-Energy Solution: Boundary pressure 
is that for which the conserved top-hat energy 
is the minimum possible for an isothermal sphere 
of fixed mass within a finite truncation radius. 

2) The Self-Similar Spherical Cosmological Infall 
Solution (Bertschinger 1985) confirms this choice 
if we identify the virialized object with the spher- 
ical region of post-shock gas and shell-crossing 
dark matter =^ explains dynamical origin of bound- 
ary pressure adopted above as the result of ther- 
malizing the energy of infall. 



Summary of the TIS Solution 



• Top-hat perturbation =^ unique, nonsingular TIS 

(minimum-energy configuration) 

=^ universal, self-similar density profile for the 
postcollapse equilibrium of cosmic structure 

— Unique scale and amplitude set by top-hat mass 
and collapse epoch 

— Same density profile for gas and dark matter (no 
bias) 

I. Matter-Dominated Case (see Table 1 and Fig. 1) 

— Finite core size: ro = 0.034 x radius n 

— Central density: po = 514 x surface density pt 

T = 2.16 Tuniform spliere — 0.72 Tsjngular isotliermal spliere 

— At intermediate radii, p drops faster than 

II. Flat, A 7^ Case 

— Profile varies with epoch of collapse, approach- 
ing case I above for early collapse. 

For example: for Qq = 1 - Ao = 0.3 for Zcoii = 
(0;0.5; 1): 

* n/ro = (30.04; 29.68; 29.54) 

* po/pt = (529.9; 520.8; 517.2) 

* T/Tuniform sphere = (2.188; 2.170; 2.163) 



Table 1: The Postcollapse Virial Equilibrium Resulting 
from the Collapse of Top-Hat Density Perturbations: 

Einstein-de Sitter Universe 



Uniform Singular Our 

Sphere Isothermal Solution* 

Sphere 

7] = ^ 05 0.417 0.554 

^M^... 1 3 2.16 

/ 2 GMm \ 

\ 5 Tvir / 

^ 1 oo 514 

Pt 

^ 1 3 3.73 

pt 

— - NA - oo 29.4 

IStt^ 18^^ (^\\^^ 

PbitcoW) 

^ 178 ^ 307 130.5 



* Our solution = minimum-energy, truncated, nonsin- 
gular, isothermal sphere 



Note: pb = cosmic mean matter density 



Density Profile of Halo which Forms from Top-Hat 



Perturbation Collapse and Virialization 
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Fig.l: Density profile of truncated isothermal sphere 
which forms from the virialization of a top-hat density 
perturbation in a matter-dominated universe. Radius r 
is in units of rm - the top-hat radius at maximum expan- 
sion, while density p is in terms of the density psus of the 
standard uniform sphere approximation for the virialized, 
post-collapse top-hat. (TIS = our solution, SUS = uni- 
form sphere, SIS = singular isothermal sphere). Bottom 
panel shows logarithmic slope of density profile. 



Direct Comparison with NFW Profile 
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Fig. 2: Continuous line = TIS profile; Dashed lines 
= "NFW" = Navarro, Frenk, and White (1996, 1997) 
profile: 



p(r) = 



5cphQ 



X 



cx{cx + 1)2' r2oo 

Range of c appropriate for X-ray clusters to early forming 
dwarf galaxies. 



Dwarf Galaxy Rotation Curves 



Q: How well does our TIS profile match the observed 
mass profiles of dark-matter-dominated dwarf galaxies? 
The observed rotation curves of dwarf galaxies can be 
fit according to the following density profile with a finite 
density core (Burkert 1995): 

PO^Burkert 



p(r) = 



(r/rc+l)(r2/r2-hl) 



A: The TIS profile gives a nearly perfect fit to the Burk- 
ert profile, (see Fig. 3) 
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Fig. 3: Rotation Curve Fit. Best fit parameters: 
Po^E!^^ = 1.216, = 3.134 

PO,TIS rO,TIS 

Solid line = Best fit TIS; Dashed line = Burkert profile, 
where a^j^ = {v^)/3 = kBT/m. 

Q: How well does this best fit TIS profile predict the 

fmax and Vmax 

f'max, Burkert ^ '^max, Burkert 1 Ql 



rmax,TIS '^max.TIS 

(i.e. excellent agreement) 



The Vmax — rmax relatioD for dwarf and LSB galaxies. 



Q: Can the TIS halo model explain the observed correla- 
tion of Vmax and Vmax for dwarf spiral and LSB galaxies? 

A: Yes, when the TIS halo model is combined with the 
Press-Schechter model which predicts the typical col- 
lapse epoch for objects of a given mass (i.e. the mass 
of the Icr-fluctuations vs. Zcoii)- (See Fig. 4) For the 
three untilted CDM models plotted, a cluster normal- 
ized Einstein-de Sitter model, and COBE-normalized 
low-density models (Qo = 0.3 and Ao = or 0.7), only 
the flat models yield a reasonable agreement with the 
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Fig. 4: Dwarf galaxies (triangles) and LSB galaxies 
(squares) from Kravtsov et al. (1998); Burkert: fit to 
data (Mori &i Burkert 2000); SCDM: Qo = 1, Ao = 
0, (TQh i = 0.5 (cluster normalized); OCDM: Qo = 0.3, 
Ao = (COBE normalized); ACDM: Qo = 0.3, Ao = 0.7, 
(COBE normalized); h = 0.65 for all. 



Galaxy Halo M — Relation 



Q: How well does our TIS halo model scaling relation 
predict the velocity dispersion of galactic halos which 
form in the CDM model according to N-body simula- 
tions? 

A: Antonuccio-Delogu, Becciani, &i Pagliaro (1999) used 
an N-body treecode at high-res (256^ particles) to simu- 
late galactic halos in the region of a single and a double 
cluster. The agreement with the TIS model is good. 
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Fig. 6: Velocity dispersion vs. mass for galactic haloes 
in cluster regions: (upper panel) double cluster, (lower 
panel) single cluster. 



X-Ray Cluster Scaling Relations 



Q: How well does the TIS halo model predict the internal 
structure of X-ray clusters found by gas-dynamical/N- 
body simulations of X-ray cluster formation in the CDM 
model? 

A: As shown below and in Fig. 5, our TIS model predic- 
tions agree astonishingly well with the mass-temperature 
and the radius-temperature virial relations and integrated 
mass profiles derived from numerical simulations by Evrard, 
Metzler and Navarro (1996). Apparently, these simula- 
tion results are not sensitive to the discrepancy between 
our prediction of finite density core and the N-body pre- 
dictions of a density cusp for clusters in CDM. 

• Mass Profile - Temperature Relation 




X/IOOO 

Fig. 5 (triangles) fit to CDM simulation results by 
Evrard, Metzler and Navarro (1996); (continuous 
line) TIS prediction. 



Mass-Temperature and Radius-Temperature Virial 
Relations 

- Evrard, Metzler and Navarro (1996) 
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- Our solution 
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/3-fits to X-ray Cluster Brightness and 
Density Profiles 
po J- ^0 

"'^^ ~ (i + rVr?)'"^'' ~ (1 + eye^^f"-'^^ 

Q: How well does the TIS model for the internal struc- 
ture of X-ray clusters predict the observed and simulated 
X-ray brightness profile of clusters? 

A: It predicts gas density profiles and brightness profiles 
which are well-fit by a /3— profile, with /3-values for the 
TIS /3-fit which are close to those of simulated clusters 
in the CDM model, but somewhat larger than the con- 
ventional observational result that f3 ^ 2/3. However, 
recent X-ray results suggest that the true /3-values are 
larger than 2/3 when measurements at larger radii are 
used and when central cooling flows are excluded from 
the fit. 

Brightness profile observations p 

Jones and Foreman (1999) 0.4-0.8, ave. 0.6 

Jones and Foreman (1992) ~ 2/3 



Balland and Blanchard (1997) 


0.57 (Perseus) 




0.75 (Coma) 


Durret et al. (2000) 


0.53 (incl. cooling flow) 




0.82 (excl. cooling flow) 


Vikhlinin, Forman, &i Jones (1999) 


0.7-0.8 


(fit by Henry 2000) 




TIS /3-fit {rc/ro,Tis = 2.639) 


0.904 


Gas density profile simulations 




Metzler and Evrard (1997) 


0.826 (DM) 


0.870 (gas) 


Eke, Navarro, and Frenk (1998) 


0.82 


Lewis et al. (1999) (adiabatic) 


~ 1 


Takizawa and Mineshige (1998) 


~ 0.9 


Navarro, Frenk, and White (1995) 


0.8 


TIS fj-f\t {vc/ro^Tis = 2.416) 


0.846 



X-ray Cluster Gas Entropy 



Q: Can the TIS model for the internal structure of X- 
ray clusters explain the observed correlation between the 
gas entropy near the cluster center and cluster virial 
temperature? 

A: Yes, but only for high-T clusters (i.e. T > few keV) 
for which energy release feedback effects were probably 
not big enough to alter the entropy of the equilibrium 
halo. (See Fig. 7) 

S = T/np, at r = O.lryir 
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Fig. 7: Cluster entropy vs. temperature. Data: Pon- 
man, Cannon, and Navarro (1999) Nature, 397; Error 
bars: T - 90% confidence level, entropy - span of varia- 
tion from r = O.Obryir to r = 0.2ryir. Our solution: thick 
line = S at r = O.lryir, dashed lines = S at r = 0.05r^;ir 
(lower), and 0.2ryir (upper). 



Cluster Mass Profiles Deduced from Strong Gravitational Lensing 



Q: Can the TIS halo model explain the mass profile 
with a finite density core measured by Tyson, Kochan- 
ski, and Dell'Antonio (1998) for cluster CL 0024+1654 
at z = 0.39 using the strong gravitational lensing of 
background galaxies by the cluster to infer the cluster 
mass distribution? 

A: Yes, the TIS model not only provides a good fit to the 
shape of the projected surface mass density distribution 
of this cluster (see Fig. 8) within the arcs, but when 
we match the central value as well as the shape, our 
model predicts the overall mass, and a cluster velocity 
dispersion in close agreement with the value ay = 1200 
km/s measured by Dressier and Gunn (1992). 
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CL 0024+1654 (Tyson et al. 1999) 
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Fig. 8: Projected surface density of cluster CL 0024+1654 
inferred from lensing measurements, together with the 
best fit TIS model. 



Summary 



• TIS profile fits dwarf galaxy rotation curves; com- 
bined with Press-Schechter formalism matches re- 
sults for observed Vmax-rmax relation for dwarf galax- 
ies 

• Predicted mass-velocity dispersion relation agrees 
with high resolution N-body simulations of galactic 
halo formation by Antonuccio-Delogu et al. (1999) 

• Predicted mass-radius-temperature scaling relations 
match simulation results from X-ray clusters in CDM 
model 

— Our solution derives empirical fitting formulae 
of Evrard, IVietzler and Navarro (1996) 

— Agrees well with X-ray cluster observations at 

• Fits high temperature X-ray cluster entropy data 

• X-ray brightness profile is predicted to match /9-fit 
with [3 ^ 0.9, larger than typically observed, but 
similar to results of gas-dynamical/N-body simula- 
tions of X-ray clusters in CDM model 

• Fits the cluster mass profile with finite core derived 
from strong gravitational lensing data of Tyson et 
al. (1999) on CL 0024-hl654 

• Predicted mass profile is close to NFW profile for 
low values of concentration parameter, outside the 
core 

This work was supported in part by NASA Grants NAG5-2785, 
NAG5-7363, and NAG5-7821, NSF Grant ASC-9504046, and Texas 
Advanced Researcin Program Grant 3658-0624-1999, and benefit- 
ted from PRS' participation at the Aspen Center for Physics in 
summer 1998 and 1999. 



References 

Antonuccio-Delogu, Becciani, &i Pagliaro (1999) To appear in the 
proceedings of tine IVIarseille IGRAP99 conference "Clustering at 
Higli Redsliift", |astro-pli/991U2bb| 

Balland C, and Blanclnard A. (1997), ApJ, 487, 33 

Bertsclninger E. 1985, ApJS, 58, 39 

Burkert A. 1995, ApJ, 447, L25 

Dressier and Gunn 1992, ApJS, 78, 1 

Durret et. al. 2000, |astro-pn/0002bi9 



Eke V.R., Navarro J.F., and Frenk C.S. 1998, ApJ 503, 569 
Evrard A.E., Metzler C.A., and Navarro J.F. 1996, ApJ, 469, 494 
Henry P. 2000, |astro-ph/UUU23bb| 



Iliev Shapiro 2000, in preparation 

Jones C, and Forman W., 1992, in Clusters and Superclusters of 
Galaxies, ed. A. C. Fabian Kluwer Academic Publishers, pp 49-70 

Jones C, and Forman W. 1999, ApJ, 511, 65 

Kravtsov A. V., Klypin A. A., Bullock J. S., and Primack J. R. 
1998, ApJ, 502, 48 

Lewis, G. F., Babul A., Katz N., Quinn T., Hernquist L., Weinberg 
D.H., |astro-ph/yyo/09/ 



Metzler C.A., and Evrard A.E., |astro-ph/9/ 10324 



Mori M., &i Burkert A. 2000, |astro-ph/0001422 



Navarro J.F., Frenk C.S., and White S.D.M. 1995, MNRAS, 275, 
720 

Navarro J.F., Frenk C.S., and White S.D.M. 1996, ApJ, 462, 563 

Navarro J.F., Frenk C.S., and White S.D.M. 1997, ApJ, 490, 493 

Ponman T.J., Cannon D.B., and Navarro J.F. 1999, Nature, 397, 
135 

Shapiro, Iliev, and Raga 1999, MNRAS, 307, 203 

Takizawa M., Mineshige S. 1998, ApJ, 499, 82 

Tyson J. A., Kochanski, G.P., and Dell'Antonio LP. 1998, ApJ, 498, 
L107 



Vikhlinin A., Forman W., &i Jones C. 1999, |astro-ph/990b200 



